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I. INTRODUCTION 



The understanding of hot dense nuclear matter is an interesting problem in theoretical 
physics in the context of heavy ion collision experiments as well as big bang cosmology. The 
problem here is basically nonperturbative. The interaction of nucleons which may arise as a 
residual interaction due to their substructure of quarks and gluons is technically not solvable. 
The usual approach is to tackle the problem through meson nucleon interactions. This 
also entails a nontrivial technical problem with the pion nucleon coupling G%j^^/4:t: being 
as large as 14.6 making any perturbative calculation unreliable. An alternative consistent 
theoretical frame work has been developed by Walecka^ consisting of interactions of nucleons 
with a neutral scalar field a as well as a neutral vector meson cu. This has been done at 
zero temperature as well as finite temperatures.^ Variations of the same model has also 
been considered including cubic and quartic terms in the a fields to reproduce correct bulk 
modulus of nuclear matter.^ These calculations however use meson fields as classical, and, 
use a (j-field which is not observed. 

An alternative model for infinite nuclear matter consisting of interacting nucleons and 
pions was considered in Ref [4]. The nuclear matter was dressed with off mass shell pion 
quanta. The scalar isoscalar pion condensates simulated the effects of a mesons^ with the 
short distance repulsion arising from composite structure of nucleons and/or through vector 
meson exchanges. This is not only aesthetically appealing with clasical a fields arising 
from quantum mechanical structures but also has a stronger phenomenological appeal as a 
mesons have not been found in nature. With a similar approach we shall consider nuclear 
matter at finite temperatures. The methods of thermofield dynamics^ fit naturally for this 
purpose because here statistical average is done through an expectation value over a "thermal 
vacuum" in an extended Hilbert space. 

The outline of the paper is as follows. In Sec. II we shall consider nuclear matter 
with pion condensates at zero temperature. Here we shall briefly summarize the results of 
Ref. [4] and show that these results constitute a particular approximation of the present 
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method. In Sec. Ill we shall discuss finite temperature nuclear matter using the methods of 
thermofield dynamics.^ We shall calculate thermodynamic quantities like pressure, entropy 
density and energy density for nuclear matter and derive the equation of state for the same. 
Section IV will consist of discussions of the results obtained in the present model. A possible 
experimental signature resulting from the model is discussed in Sec. V. In Appendix A we 
shall summarise some results of thermofield dynamics^ for the sake of completeness. 

II. ZERO TEMPERATURE FORMALISM 

We shall here start with the effective Hamiltonian for pion nucleon interactions which 
was derived in Ref. [4] and given as 

nN{x) = nP{x) + Hint{x) (1) 

where the free nucleon part H^^ is given as 

nP{x)=,p\{x)e,M^) (2) 

and the effective pion nucleon interaction part is given as 

iG^^, G^ 



M^)- (3) 



In the above, — (M^ — S/IY^"^ where M denotes the nucleon mass. Furthermore the free 
meson part of the Hamiltonian is given as 

^m(^) = ^ [4>'l + {\/<pi) ■ (V0i) + m^<i>f\ (4) 

Clearly in the above i^ij refers to the large component of the nucleon spinor and = rj0j. 

We expand the field operator (pi^z) in terms of the creation and annihilation operators 
of off-mass shell mesons satisfying equal time algebra as 

(t>i{z) = -^={ai{z)^ + ai{z)) (5) 

and 



M^) = - o^^)) (6) 

where we take with the perturbative basis Uz — {w? — Vz)^^^) with m denoting the mass of 
the meson. 

The two pions in Eq.(3) constitute a scalar-isoscalar interaction of nucleons and thus 
could simulate the effects of cr-mesons. With this in mind, let us consider a two pion 
creation operator given as 

- \ j f{k)ai{k)^a,{-k)Uk (7) 

where the arbitrary function f{k) will be determined by a variational procedure. We now 
introduce a 'meson' dressing of nuclear matter through the state 

1/ >= U\vac > (8) 

where 

U = exp{B^ - B) (9) 
In contrast to Ref . [4], the operator U is unitary, and thus we have 

</!/>=! (10) 

This makes the formalism simpler and more realistic. Further with this ansatz, one easily 
obtains 

U^ai(k)U = (coshf(k))ai{k) + {sinhf{k))ai{-ky (11) 

which is a Bogoliubov transformation to be used later for the calculations. 

We shall take nucleons occupying a spherical volume V — 47r^, where N/{A'k^) — p 
remains constant as iV — > oo and we neglect the surface effects. We describe the system 
with the density operator pjv such that^ 

^r[pjv'0/3(y) Va(^)] = Pa/six, y) (12) 
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and 

tr[pNN] = J paa{x, x)dx = N = pV (13) 
We obtain the total nucleon energy as 

hf = <f\tr[p^nP{x)]\f> 



7 , 3 k 



-—ki(M + ^) (14) 

67r2 lOM^ ^ ^ 

where 7 = 4 for nuclear matter and p and kf are related by the equation p = 7A;^/67r^. With 
the meson field operator expansion as in Eq. (5) and (6) we may write eq. (4) as 

Hm{x) = ai{xyLj^ai{x) (15) 

Using Eq.(ll) we now obtain kinetic energy density due to the mesons as 

hk = < f\nM{x)\f > 

= 3(271)-^ J dku;{k)sinh^f{k) (16) 



where uj{k) = y k"^ + rn?. We next proceed to evaluate from Eq.(3) the interaction energy 
density, with ex — M, 

hint = < fltrlpNT-Cintix)]]/ > 

^% < /I : 0i(x)0,(f) : 1/ > (17) 



2M 

Using the Bogoliubov transformation (11) we have from equation (|l 

G'^p ,_o f dk I sinh2f{k) . „ 



hut = ^ ■ 3(2vr)-^ j ( ^ + s^nh'f{k) ) (18) 

The meson energy density thus is given as 

hm = hint + hk (19) 
Now extremising equation (|TUp with respect to f{k) we obtain the solution 
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tanh2m = ■ (20) 



2M uj%k) + 



2M 



We may compare the same with the results of Ref. [4] which corresponds to the first term 
in the expansion of left hand side of Eq. (120|). We then obtain the corresponding meson 
energy density from the kinetic and interaction terms as 

hm = hk + hint 

3 .-s.G^p^ r dk 



2 ^ ^ ^2M' J 



uj{k)[u{k){u{kY + ^y/^ + {ujikf + )] 

(21) 



We note that Eq.(21) is not acceptable since the energy density diverges. This happens 
because we have taken the pions to be point like and assumed that they can approach as 
near each other as they like, which is physically not correct. If we bring two pions close to 
each other there will be an effective force of repulsion because of their composite structure. 
We thus assume a phenomenological term corresponding to meson repulsion as 

h^ = 3a(27r)-3 j {sinh^ f{k))e^'''' dk (22) 

where a and R.,, are two parameters to be determined later. Now extremising equation (16), 
(18) and (22) with respect to f{k), we obtain 

tanh2f(k) = ^ ^ -5 ^ . (23) 

2M uj^(k) + ^ + auj{k)e^i''' ^ ' 

In place of Eq. (21) we now obtain the expression for hm as 



hm — hk + hint + h. 



3 ^-3 ^G'P 



(2 



TT 



dk 



2 ^ ' '2M' 



1 



(24) 

where lo = u{k). Finally we have to include the energy of repulsion which may arise from 
vector meson interaction and/or from finite size of the nuclei. We shall here parametrize the 
effect of such a repulsion contribution by the simple form 



Hr = Ap" (25) 

where A is an arbitrary constant to be fixed from phenomenology. We note that equation 
(25) can arise from a Hamiltonian density given as 

'Hij(f) ='0(f)V(^) j vr{x - y)ip{y)'^ip{y)dy (26) 

where, when density is constant, we infact have 

A = y VR{f)dr 

We next minimise the energy per nucleon given by 

^ ^ {hm + hf + Hr) ^27) 
P 

as a function of p, with the parameters a, and A to be subsequently fixed. Extremizing 
the single particle energy, for A = 0.54 fm^, = 1.18 fm and a = 0.12 GeV we have 
the single particle energy given as £^ = —15.03 MeV at the saturation density p — 0.19 
fm^^corresponding to kf — 1.42 fm~^. The equation for E vs. p is shown in Fig. 1 . The 
incompressibility of the nuclear matter is given as 

K = k)—j = 151.2MeV (28) 
dkf 

We may note that the parameters and the results of the present analysis do not differ 
significantly from those of Ref. [4]. As may be seen in Eq.(20) or (23) Ref. [4] is an 
approximation of the present framework. 

We may also calculate the average pion number per nucleon p^ given as 

Pm = < /|al(^)ai(^)|/ > 

= 3(27r)-^ J dk sinh^f{k) (29) 

which we plot in Fig. 4 as a function of nucleon density. This may be relevant for heavy ion 
collisions where nuclear matter gets compressed. 
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III. FINITE TEMPERATURE FORMALISM 



We shall now generalize the formahsm developed in Sec. II for finite temperatures. For 
this purpose we shall use the methodology of thermofield dynamics.^ In Appendix A, we 
summarize the salient features of the same. Here Eq.(7) for the background pion will become 
modified with introduction of extra thermal modes as illustrated in the Appendix. In fact, 
we shall have the temperature dependent background off-shell pion pair configuration given 
as 

>= Uimf >= Ui{P)U\vac > (30) 

where U is the same as in Eq.(9), and ^7/(/3) describes the change with temperature. The 
expression for this in terms of ordinary and thermal modes is given as 

Uj((3) = exp(B\((3) - Bj{(3)) (31) 

where, as in Eq.(A4) of Appendix A, 

S,(/3)t = 1| eB{k,(3Mk)%{-k)Uk (32) 

In the above, 

bi(ky = Uai(k)U^ (33) 

— * ~* J. 

so that bi{k)\f >— 0. Thus bi{ky creates excitations over the zero temperature configuration 
given by 1/ >. As before we shall calculate the energy expectation values, but here shall 
minimise the thermodynamic potential as appropriate at finite temperatures. To do so we 
note that a little algebra yields 

U^ai{k)U = (coshf(k))(cosheB(k,(3))bi{k,(3) + (sinhf(k))(cosheB(k,(3))bi{-k,(3y 

+ {sinhf{k)){sinheB{k,P))bi{-k,P) + {coshf{k)){sinheB{k,P))bi{k,py 

(34) 
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which again is a Bogohubov transformation, and will be used for obtaining the energy 
expectation values. The parallel unitary transformation as in Eq.(30) for the temperature 
dependance in nucleon sector for fermions is given as 

UniP) = exp{B\j{(3) - Bjj{P)) (35) 

with 

Bnipy = i I Opik, l3)tl^i{k)^^i{-k)Uk (36) 

where 0F{k,j3) will be determined later. As earlier for thermal averages, we shall replace 
PnN in Eq. (13) by 

PNMx)Ha{x) ^ Un{P)Ha{x)Ha{x)Unm. 

We then have the nuclear matter density 

p{(5) = < vac\Uu{P)Ha{x)H^{x)Uum'^(^c > 

= 'y{2n)-^ J dhsin"^ Of (37) 

As noted in the Appendix A, sin^Op is the distribution function for the fermions. Clearly, 
with sin^Op — Q{kf — k) Eq. (37) gives p — ^kj/Gn'^ of zero temperature. sin^Op for 
the interacting system will be determined here from the construction of the thermodynamic 
potential. 

We shall now calculate different contributions to the energy expectation values corre- 
sponding to the Hamiltonian as in Eq. (1) and (4) as well as Eq. (25). We thus have for 
the nucleon kinetic term 

hf{l3) = < vac\Uim^i,i{x)^^^^^i,i{x)UHmvac > 

^^(27,)-' j dk^sin-'ep (38) 

The kinetic energy due to the mesons is given by 
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hk{p) = <f,p\nM{x)\f,P> 

= 3{27r)-^ J dkuj(k)[sinh'^f(k)cosh2eB(k,l3) + sinh'^OBik, P)] (39) 



where uj{k) = y k'^ + m?. In the above, we may note that when ^ it reduces to the 
Eq. (16). We next derive the interaction energy density from Eq. (3) as 

h.UP) = < f,l3\UimHi{x)Hi{x)^cl>\x)Un{l3)\f,l3> 

2tx 
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^p(/3)</,/?|:0,(f)0,(f) :!/,/?> 



(40) 



where 



3 f dk sinh2f(k)cosh26B . ,2?/ A r ■ iin \ /.-,\ 

hM = jTT^ / — ^ + sinh^f{k)cosh2eB + smh^OB • (41) 

{2Txy J uj{k) \ 2 I 



As in Eq. ( p2D we shall now assume a phenomenological term corresponding to meson 
repulsion due to composite structure of mesons given as 

h^{(5) = 3a{2'Ky^ J (sinh^f{k)cosh2eB + sinh^ee) e^'^^'dk (42) 

Finally, the nucleon repulsion term parallel to Eq.(25) is 

hR = Xp\P) (43) 

where p(/3) is as given in Eq. (37). Thus the energy density is given by 

E{(3) = {hf{f3) + hUf^) + hn{f3))/p{(3) (44) 

where 
as before. 

The thermodynamic potential density Q is given by 

n{f3) = E{(3)p - ^ - /iP (45) 
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where the last term corresponds to nucleon number conservation with /i as the chemical 
potential. We may note that we shall be considering temperatures much below the nucleon 
mass so that in the expression for we do not include antiparticle channel. The entropy 
density above is a — ap + ctb with ap being the entropy in fermion sector given as^ 

ap = -7-^ / dk\sin^eFik,l3)lnisin'^eFik,l3)) + cos^eF(k,l3)lnico^ . 
and simillarly the meson sector contribution ctb is given as 

— "TTTT^ ( dk sinh?9B{k-i P)ln{sinh?9B{k^ (5)) — cosh^9B{k, P)ln{cosh^9B{k, (5)) . 
Thus the thermodynamic potential density now is a functional of 9F{k,(3), 9B{k,l3) as well 

~ — * 

as the pion dressing function f{k) which will of course depend upon temperature. Extrem- 
isation of Eq (45) with respect to f{k) yields 

tanh2f(k) = ^ ^ (46) 

which is of the same form as Eq. (23) for zero temperature case except that now p is tem- 
perature dependent through Eq. (37). Similarly minimising the thermodynamic potential 
with respect to 9B{k,P) we get 

sinh'9B = ^^J—^ (47) 

where 

co' = (uj+ 4^ + ae'^'^')cosh2f(k) + -^sinh2f(k) (48) 
Once we substitute the optimised dressing as in Eq. (46), the above simplifies to 

^' ^ (^ + ^ + ae^^^y/^^^ ae'^'^'y^ (49) 

which is different from u) due to interactions. Further, minimising the thermodynamic 
potential with respect to 9F{k,P) we have the solution 



with 

e. = ^/2M + 2Ap+^ (51) 

where I2M is given in equation (41). We may note that the change in ep above from k'^ /2M 
is also due to interaction. 



IV. RESULTS 

To calculate different thermodynamic quantities as functions of baryon number density 
we first use Eq. (37) to calculate the chemical potential in a self consistent manner with 
p and 11 occuring also inside the integrals through sin^Op as in Eq. (50). Thus for each 
p, we determine p so that Eq. (37) is satisfied. The ansatz functions 9b, Op and f{k) get 
determined analytically through the extremisation of the thermodynamic potential and the 
parameters a, A and Rt^ are as obtained in Sec. II. 

With the thermodynamic potential determined as above, we calculate different thermo- 
dynamic quantities. We first plot the binding energy per nucleon as a function of baryon 
number density in Fig. 1 for temperatures zero MeV, 5 MeV, 10 MeV and 15 MeV. As 
the temperature increases the minimum shifts towards higher densities. This may be un- 
derstood from the fact that to compensate for the larger kinetic energy a large value of p 
is needed to give the minimum until the nuclear repulsion effects take over and hence again 
energy increases. We may compare our curves with those of Ref [2] . Our curves compared 
to theirs are rather shallow. This is consistent with the fact that we get a smaller value of 
compressibility (~ 150 MeV). 

We next plot pressure as may be defined from thermodynamics^ 

P((3) = (52) 

The same is plotted in Fig. 2. Vanishing of the pocket in the above curves for temperatures 
around 15 MeV may be noted which could be an effect of phase transition. 
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Entropy density is plotted as a function of nucleon density in Fig. 3 at temperatures 50 
MeV, 100 MeV, 150 MeV and 200 MeV. The general tendency appears to be quite similar 
to that of Ref. [2]. 

Finally, in Fig. 4 we plot the average number of pions per nucleon (p^) as a function 
of nucleon density p. It is interesting to note that this result is practically independent 
of temperature. Therefore only the zero temperature graph is plotted in Fig. 4. For low 
nucleon density pm rises rather sharply and almost levels off at p > 0.3 fm""^ ~ x 1.5 
where ps is the density of stable nuclear matter at zero temperature. We note that the rise 
in meson number with compression may be observable in heavy ion collisions. 

V. DISCUSSIONS 

In the present model we have considered nuclear matter at zero temperature and finite 
temperatures in a nonperturbative manner for the interacting pion nucleon system using 
thermofield method. This yields results similar to those of mean field approximation, with 
(7-mesons not being needed, and, the calculations being quantum mechanical. 

Besides the aesthetic appeal, we note that it may have experimental consequences. For 
example, the off shell tt'^tt" pair as present through pion dressing may annihilate to hard 
photons with a probability in excess of what one may expect otherwise. This is likely to 
depend on individual nuclei which we may dress with pions as here. 
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APPENDIX A 



We here summarize briefly the sahent features of thermofleld dynamics^ as used in the 
present paper. 

In statistical mechanics, the thermal average of an operator O is given as, with (5 = 1/kT, 

< O (Al) 

where, the trace is taken over a complete basis of states. First we note that in the zero 
temperature limit, the above reduces to ground state expectation value for the operator O. 
This is easily seen as 

A < I 6 I > e-^'°+ < 1 I 6 I 1 > e-^^^ + ■■■ 
lim < O >g — hm ^ s 

/3-*oo ^ /3-»oo e-P^o + e"^^i H 

<0|0|0> + <1|6|1> e-^(^i-^°) + • • • 

= lim ^7 r 

^^oo 1 + e~'^('^i~^o) H 

= < I O I >, (A2) 

where | > corresponds to the state with the lowest energy. In thermofleld method, one 
essentially generalises (742) to the case of flnite temperature and deflnes a "thermal vacuum" 
such that the statistical average reduces to an expectation value with respect to the thermal 
vacuum. Thus we want that for some | 0(/3) > the relationship 

< O >/3= j.\^-f,H^^ ^< m I O I 0(/?) >, (A3) 

where | 0{j3) > is deflned as the "thermal vacuum". This can be done if one doubles the 
degrees of freedom i.e. corresponding to every physical operator a, a"tilde" operator d is 
introduced. In the above, d{k) are the new operators named as "thermal modes". They 
are associated with negative energy, with conventional quantisation, but, do not have any 
physical signiflcance in the sense of observation of these modes. In zero temperature vacuum, 
these modes are absent, so that conventional fleld theory holds. At flnite temperature the 
ground state is replaced by | 0(/3) > given as 
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|0(/3) > =UB{P)\vac> 

= exp{ j eB{k,l3){a{k)^a{-k)^ - h.c.)dk) \ vac >, (A4) 

— * , 

where a{—ky in the above corresponds to the extra Hilbert space. It is now convenient to 
define a thermal basis 



^ a{k,/3) ^ 



a{-k, 



U{P) 



( a{k) ^ 



a{-ky 



U{l3)-\ 



(A5) 



which amounts to the Bogohubov transformation 



a{k,(5) 
d(-k,py 



^ ( cosheB{k,(3) -smheB{k,(3)^ ^ 



V 



a{k) 
a{-ky 



(A6) 



— sinh 6b {k, (5) cosh 6b {k, (5) 

a{k,(3) and d^{k,P) are the annihilation and creation operators at temperature (3 — 
corresponding to the thermal vacuum such that 

a{k,P) I 0(/3) >= = a\k,P)W) > . 

The function 6B{k,(3) is calculated through minimization of thermodynamic potential den- 
sity given as 



(A7) 



where /i, the chemical potential corresponds to a conserved number, and the entropy density 

3 



(Jb = 



dk 



sinh^6B{k, p)ln{sinh^6Bik, f3)) - cosh'6B{k, P)ln{cosh'6B{k, f3)) 



For zero chemical potential and for free fields, extremization of the free energy then yields 

1 



sinh 6b — 



where we have taken the Hamiltonian density as = / {z)ujza{z)dz so that 

E((3) = < vac, (3\n°\vac(3 > 
1 



(A8) 



sinh 6Bi^{k)dk 
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If we substitute this value the free energy density becomes 

which is the same as derived through a more convenational treatment of temperature de- 
pendent quantum field theory.^ As stated, we do here the analysis through thermofield 
dynamics since the calculations are simpler. In case of interacting field however the solution 
for 9B{k,P) will not be given by equation (A8) and will depend upon the interaction. 
Similarly for the fermionic sector the thermal vacuum will be given by 

|0(/3) > = UF{P)\vac> 

= exp( J epik, (3) ('ijji(k)^i(-k)'' - h.c.)dk) \ vac > (A9) 

where ^\ corresponds to the creation of the fermionic thermal modes. The entropy is given 



as 



Up = 
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j dk[sin^eF{k, (3)ln{sin'^eF{k, (3)) + cos^epik, (3)ln{cos'^eF{k, (3)) 



The thermal Bogoliubov transformation is now given as 



cos^f(A;,/3) —sm9p(k, P) 
— sin 9F{k, (3) cosOpik, 13) 



J 



(AlO) 



Parallel to equation (A8) the function of 9F{k,(3) is given as 

1 



(All) 



where n is the chemical potential corresponding to baryon number conservation and 



uj{k, P) — \/W~^rM^ for free fermions of mass M, so that when statistics is known, the 
corresponding Bogoliubov tranformation relating zero temperature ground state with the 
thermal ground state of the extended Hilbert space is known. The ground state or the 
"thermal vacuum" obviously contains particles with appropriate distributions. 

The methodology enables us to replace mixed states of statistical mechanics by pure 
states in an extended Hilbert space while generating correct distribution functions. The 
extra thermal modes enable us to do this. 
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FIGURES 

FIG. 1. Binding energy per nucleon E as a function of nuclear matter density p at temperatures 
T=0 MeV, 5 McV, 10 MeV and 15 MeV. Tlie gradual shift in the energy minimum towards higher 
densities may be noted. 

FIG. 2. Pressure P as a function of nuclear matter density p at temperatures T=0 MeV, 5 
MeV, 10 MeV, 15 MeV and 20 MeV. The vanishing of 'pocket' at about 15 MeV may be taken as 
a signature of phase transition. 

FIG. 3. Entropy density cr as a function of nuclear matter density p at temperatures T=50 
MeV, 100 MeV, 150 MeV and 200 MeV. 

FIG. 4. Average number of pions per nucleon p^ as a function of nuclear matter density p. 
Since the results are effectively independent of temperature, only the zero temperature graph is 
plotted. 
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